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Like many others, silicate solids dissolve when placed in contact with water. In a given aqueous 
environment, the dissolution rate depends highly on the composition and the structure of the solid, 
and can span several orders of magnitude. Although the kinetics of dissolution depends on the 
complexities of both the dissolving solid and the solvent, a clear understanding of which critical 
structural descriptors of the solid control its dissolution rate is lacking. Through pioneering disso¬ 
lution experiments and atomistic simulations, we correlate the dissolution rates - ranging over four 
orders of magnitude - of a selection of silicate glasses and crystals, to the number of chemical topo¬ 
logical constraints acting between the atoms of the dissolving solid. The number of such constraints 
serves as an indicator of the effective activation energy , which arises from steric effects, and prevents 
the network from reorganizing locally to accommodate intermediate units forming over the course 
of the dissolution. 


I. INTRODUCTION 

Due to their abundance in Earth’s crust, silicate min¬ 
erals and glasses are the base components of numerous 
industrial materials including glasses, ceramics, and con¬ 
crete. As such, understanding and predicting the kinet¬ 
ics of silicate dissolution in the presence of water has 
fundamental and industrial implications, with examples 
including bioactive glassed, nuclear waste confinement 
matrices 2 3 and the reaction of cement with water, the 
binding phase in concrete, all of which require accu¬ 
rate predictions of dissolution rates, in the context of 
microstructure evolution, phase stability and chemical 
durability. 

Predicting the kinetics of dissolution of materials under 
given thermodynamic conditions is, a priori, a complex 
problem as it depends, amongst others, on the solvent 
chemistry 5 ^ and the geometry of the surface^- 7 . Another 
difficulty arises from the fact that the dissolution of min¬ 
erals can feature leaching, or incongruency, that is, occur- 
ing in a nonstoichiometric fashion 7 . However, this latter 
feature typically is true only during the initial stages of 
contact with water, with steady-state dissolution usually 
showing congruency®! Nevertheless, for a given solvent 
and thermodynamic condition, both the composition and 
structure of a material dictate its dissolution kinetics 9 10 . 

Mere knowledge of the composition of a material is in¬ 
sufficient to predict its dissolution rate as, e.g., glassy sil¬ 
ica dissolves between one and three orders of magnitude 
faster than crystalline a-quartz^, depending on the sol¬ 
vent pH. While it may be argued that glasses should in¬ 
deed be less stable that their crystalline equivalents, such 
a trend is not systematic as, e.g., albite glass and crystal 
show comparable dissolution rated- 11 . In fact, the disso¬ 
lution of silicate glasses and minerals has recently been 


shown to occur by similar mechanisms and to be mostly 
controlled by short-range atomic order, rather than the 
long-range disorder 12 . Indeed, the dissolution rates of 
silicate glasses are generally thought to decrease with 
the connectedness, or the degree of polymerization of 
the atomic network 13 , described by the ratio of terminat¬ 
ing non-bridging oxygen (NBO) per tetrahedral unit (T) 
formed by silicon or aluminum atoms (NBO/T). How¬ 
ever, using such a metric is far too restrictive as, e.g., 
silica, albite, jadeite, and nepheline glasses all feature a 
fully polymerized network (NBO/T = 0), but show vastly 
differing dissolution kinetic^ 14 . As such, given the disor¬ 
dered nature of glass, it is a challenge to identify a met¬ 
ric that is simple enough to enable practical predictions 
of dissolution rates, but that takes into account enough 
structural detail to provide realistic results. 

By capturing the topology of atomic networks while 
filtering out chemical details that ultimately do not af¬ 
fect the_ macroscopic properties, topological constraint 
theory®^ (TCT) is a promising tool that can eluci¬ 
date such a metric. Within the TCT framework, atomic 
networks are described as mechanical trusses, in which 
the atoms experience topological constraints, as imposed 
by radial and angular chemical bonds. Hence, follow¬ 
ing Maxwell’s stability criterion^, an atomic network is 
described as flexible, stressed-rigid, or isostatic, if the 
number of topological constraints per atom (n c ) is lower, 
higher, or equal to three, i.e., the number of degrees of 
freedom per atom in three dimensions (see Figure |T|) . 
In this paper, we establish that the dissolution rates of 
glassy and crystalline silicates in caustic environments 
are dictated by n c , which serves as an indicator of the ef¬ 
fective activation energy of the dissolution process, and, 
in turn, of the activation energies of solid state ion- 
diffusion and conduction. This suggests that, simply, the 
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Figure 1. A sketch representing the three states of rigidity 
(flexible, isostatic, and stressed-rigid) of the atomic network 
of a solid in contact with water. The flexible state is the 
only one that features internal low-energy (floppy) modes of 
deformation. 



Figure 2. An illustration of the dissolution of o-quartz par¬ 
ticulates, as visualized using vertical scanning interferometry 
(VSI), before and after 90h of solvent contact at pH 13 and 
25 °C. Dissolution is tracked by measuring the decrease in 
particle height over time. 


mechanical stability of the bulk atomic network controls 
the dissolution behavior, and the transport properties of 
silicate solids. 


II. RESULTS 

A. Dissolution rates 

The dissolution rates of glassy silica and <a-quartz were 
measured under isothermal conditions, and at different 
temperatures using vertical scanning interferometry 19 
(VSI, see the Supplemental Material). This technique 
has been extensively applied to measure the dissolution 
rates of minerals of geological relevance 20 21 . By directly 
tracking the evolution of the surface topography in time, 
with sub-nanometer vertical resolution, VSI accesses the 
true dissolution rate of a solid dissolving in a given sol¬ 
vent. Unlike dissolution assessments that are based on 
analysis of solution compositions, which may be affected 
due to aspects including metastable barrier formation, 
incongruency in dissolution or, ion adsorption, VSI ana¬ 
lytics are not influenced by such complexities. Dissolu¬ 
tion rates of the solids were quantified using a rain-drop 
procedure 19 !, wherein both the solution pH and composi¬ 
tion (i.e., the under-saturation level with respect to the 
dissolving solids) are kept constant over the course of 
a given experiment (see Figure |2j. In addition to our 
measurements, to assess the compositional dependence 
of dissolution kinetics, the dissolution rates of sodium 
trisilicate (Na 2 0-3Si0 2 ), albite (Na 2 0 -Al 2 03 - 6 Si 0 2 ), 
jadeite (Na 2 0 -Al 2 03 - 4 Si 0 2 ), and nepheline (Na 2 0- 
Al 2 03 - 2 Si 0 2 ) gla sses were s ourced from the tabulations 
of Hamilton et all 11 13 14 22 . Note that, for consistent 
comparisons between solids of different compositions, all 
the dissolution rates presented herein are normalized in 
terms of moles of 0 2 dissolved per unit of surface and 
time. 


B. Relationship between composition and 
dissolution rates 

Figure^ shows dissolution rates for a range of solids in 
basic aqueous solvents. As expected, all the dissolution 
rates increase with pH, due to the increasing abundance 
of HO - species available for nucleophilic attacks of Si or 
A1 tetrahedra 22 . Although the slope of the dissolution 
rates versus pH curve varies for each material, we note 
that the ranking of dissolution kinetics remains consistent 
across a wide pH range and scales as <a-quartz < albite 
< jadeite < nepheline < glassy silica < sodium silicate. 

The structural origin of the dissolution rates can be 
elucidated by comparing dissolution rates for a represen¬ 
tative pH, in this case pH 12. First, as shown in Figure 
[ 3 ) 3 , relying on the degree of depolymerization, NBO/T, 
is insufficient. For example, if sodium silicate (NBO/T = 
0.67) dissolves faster than glassy silica (NBO/T = 0), this 
metric does not distinguish the structures of nepheline, 
jadeite, albite and <a-quartz (NBO/T = 0) either, al¬ 
though their dissolution rates span three orders of mag¬ 
nitude. Second, one cannot rely on composition either. 
The ranking of the energies of Si-0 (444 kJ/mol), Al- 
O (423 kJ/mol), and Na-0 (83 kJ/mol) bonds suggests 
that the stability of a given silicate should increase with 
the fraction of silicon present, thereby decreasing the dis¬ 
solution ratJ^ 3 . Although, as shown in Figure [ 3 ) 3 , this 
trend is true for nepheline, jadeite, and albite, it does 
not allow one to explain differences in dissolution rates 
between glassy and crystalline silica. In summary, none 
of the typical, and expected metrics are able to capture 
sufficient structural detail to permit quantitative estima¬ 
tion of the dissolution kinetics. 

We now calculate the rigidity of the atomic networks 
of the different silicate solids to evaluate the relevance 
of this metric for dissolution. In a fully connected net¬ 
work, like, e.g., chalcogenide glasses, an atom of coordi¬ 
nation number r experiences r/2 radial bond-stretching 
(BS) constraints, and 2r - 3 angular bond-bending (BB) 
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Figure 3. Influence of composition and structure on dissolution rate, (a) Room-temperature dissolution rates (K) of sodium 
silicatd^, glassy silica, nephelinJ^, jadeit^P^, albite 14 and o-quartz (normalized with respect to the amount of O 2 dissolved) 
versus solvent pH. (b) Dissolution rates at pH = 12 with respect to the ratio of non-bridging oxygen and Si or A1 tetrahedra 
(NBO/T). (c) Same data as b), plotted with respect to the ratio of Si and cations atoms (Si/(Si+Al+Na)). (d) Same data 
as b), plotted with respect to the number of constraints per atom n c . The data are fitted by an Arrhenius-like law K = 
Ko exp(— n c Eo RT ), with Kq = 1.4 10 11 /xmol/m 2 /s and Eo = 25.5 kJ/mol (dashed line). The inset shows the same data as d), 
in linear scale. The grey area indicates the flexible domain (n c < 3). 


constraints 24 . This allows one to enumerate n c = 3.67 
constraints per atom for quartz^ (stressed-rigid). In 
glassy silica, the Si-O-Si angle typically shows a much 
broader angular excursion^, so that this angle appears 
unconstrained. This results in n c = 3.00 for glassy silica 
(isostatic), in agreement with the observation that sil¬ 
ica is an excellent glass forming system 25 . The insertion 
of sodium cations typically substantially decreases the 
glass transition temperature of the glass and, thereby, 
restores the Si-O-Si BB constraint. In addition, each 
Na cation depolymerizes the base silica network by cre¬ 
ating one NBO. The number of constraints per atom of 
(Na 20 )a;(Si 02 )i-;E can then be calculated 27 ^ as being n c 
= (11 - 10x)/3. This is in agreement with experimen¬ 
tal evidence^ of a rigidity transition (n c = 3) at x = 
0.20. Then, from a simplistic viewpoint, starting from 
sodium silicate, each added aluminum cation consumes 
one NBO. Na cations then act as charge compensators 
around A1 IV cations 29 and are, therefore, not considered 
any further in the constraints enumeration. Note that the 
tetrahedral angular environment of A1 cations is not as 
well defined as that of Si, so that the 0-A1-0 angles are 
considered unconstrainecP 0 . Consequently, the number 
of constraints per atom for (Ndi 20 ) x (A\ 20 s) x (Si 02 )i- 2 x 


is given by n c = (11 - 10x)/3. However, this enumer¬ 
ation under-estimates n c as a small fraction of over¬ 
coordinated three-fold o xygen tricluster and five-fold alu¬ 
minum species are found 29 30 . Hence, the rigidity of the 
aluminosilicate glasses was evaluated by a careful analysis 
of their structure, as predicted by the molecular dynam¬ 
ics (MD) simulations. This allows clear differentiation 
of intact from thermally broken BS (BB) constraints by 
computing the radial (angular) excursion of each neigh¬ 
bor, using a previously established methodology 27 31 32 1 
(see the Supplemental Material). As such, we find n c = 
3.33, 3.22, and 3.00 for albite, jadeite, and nepheline, 
respectively. Note that the latter value implies that 
nepheline features an isostatic network (n c = 3), which is 
consistent with the fact that its dissolution rate is similar 
to that of glassy silica. 


C. Topology controls the kinetics of dissolution 

Figure[3]i shows the dissolution rates K as a function of 
the number of constraints per atom n c . As shown in the 
inset, K significantly increases in the flexible domain (n c 
< 3), in agreement with the idea that, due to the presence 
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of internal low energy modes of deformation (i.e., floppy 
modes), a flexible network is less stable than its rigid 
counterparts, and should, therefore, dissolve faster. The 
observed evolution of K with respect to the number of 
constraints n c suggests an Arrhenius-like relationship, of 
the form: 


- which we identify as the rigidity of the network. Sec¬ 
ond, the activation energies all increase with the num¬ 
ber of constraints per atom and are similar to the values 
predicted by equation [2j that is, the effective activation 
energy induced by the number of constraints per atom. 


K - Kq exp 


f n c E 0 \ 

V RT ) 


(i) 


where Kq is a rate constant that depends on the solu¬ 
tion phase chemistry, and yields the barrier-less disso¬ 
lution rate of a completely depolymerized material (i.e., 
for which n c = 0) and Eq = 25.5 kJ/mol is an energy 
barrier that needs to be overcome to break a unit atomic 
constraint. This suggests that: 


Ef = n c E 0 (2) 

serves as an effective activation energy that controls the 
kinetics of the dissolution. To validate this hypothesis, 
we evaluated the low-temperature activation energy (Ea) 
of dissolution of <a-quartz and glassy silica by measuring 
their dissolution rates at 3.5, 25, and 45 °C at pH 13. 
We chose this higher pH to achieve appreciable dissolu¬ 
tion of quartz at the lowest temperature. Note that these 
materials feature a different n c (3.67 and 3.00 for quartz 
and glassy silica, respectively) but similar composition, 
which allows us to isolate the singular effect of atomic 
topology on the activation energy, with no compositional 
effects. As illustrated in Figure [4^i, both materials show 
an Arrhenius-like evolution of their dissolution rate with 
temperature. As such, we find Ea = 101T33 kJ/mol 
and 52T20 kJ/mol for <a-quartz and glassy silica, respec¬ 
tively, in agreement with prior studies 33 34 - an expected 
outcome since <a-quartz shows a higher n c than glassy 
silica (i.e., on account of its stressed-rigid nature). Inter¬ 
estingly, as shown in Figure [4]3, these values are in fair 
agreement with the effective activation energy predicted 
simply from the knowledge of the number of constraints 
per atom (Ej^ = ti c Eq ), which yields to EJ^ = 93.6 
kJ/mol and 76.5 kJ/mol for <a-quartz and glassy silica, re¬ 
spectively. This suggests the following atomistic picture: 
starting from n c = 0, which would correspond to a fully 
depolymerized material, each new constraint per atom ef¬ 
fectively reduce the dissolution kinetics by increasing the 
associated activation energy needed for bond rupture. 

The activation energies for dissolution are compared to 
the activation energies of other thermally-activated pro¬ 
cesses, namely, for sodium silicate (self-diffusion 35 and 
conduction^ of Na), potassium silicate (self-diffusion 35 
and conduction^of K), fused silica (diffusion of wateiP^), 
and quartz (diffusion of watei®). First, we note that, in¬ 
terestingly, for a given material, the activation energies 
associated with dissolution, diffusion and conduction are 
consistent with each other, which suggests that these en¬ 
ergy barriers all arise from a common atomistic origin 


III. DISCUSSION 

The question remains regarding the origin of the re¬ 
lationship that is found herein between the activation 
energy Ea and the number of topological constraints per 
atom n c . As such, we propose that n c serves as an indica¬ 
tor of steric effects in the atomic network, which prevent 
reorganization and internal motion of the constituent 
species. The reaction mechanism of a silicate with water 
can take different forms: (1) hydration, where intact wa¬ 
ter molecules enter the glass network^, (2) hydrolysis, 
where water reacts with Si-0 bonds to form hydroxyl 
groups, and (3) ion-exchange, where network modifying 
cations, e.g. Na + , are replaced by H + cations^. The 
three mechanisms can occur simultaneously, and influ¬ 
ence process kinetics. For (1), the rate of water diffusion 
is primarily controlled by steric constraints, i.e., atomic 
packing, imposed by the network 33 . However, in typi¬ 
cal silicates, the holes in the structure, as defined by the 
ring statistics, are too small to allow direct diffusion of 
intact water molecules into the network. Therefore, wa¬ 
ter would preferentially penetrate through hydrolysis and 
condensation reactions 22 . In mechanism (2), in basic sol¬ 
vents, Si tetrahedra are subjected to nucleophilic attack 
by HO - species, thereby forming five-fold coordinated Si 
intermediate units 33 . On the contrary, in acidic environ¬ 
ments, inter-tetrahedral bridging oxygen atoms are sub¬ 
jected to electrophilic attack by protons (H + ) and form 
three-fold coordinated O intermediate units 33 . In these 
two cases, process kinetics depend on the ability of the 
network to accommodate such over-coordinated units by 
local rearrangements 33 . For (3), the rate-limiting step 
is suggested to be the penetration of the glass surface 
by watei®. However, the three mechanisms are strongly 
coupled as the depolymerization of the silicate network 
due to (2) can open up the atomic network, e.g., due to 
ring opening thus enhancing water mobility for (1) and 
(3) to occur at an expedited rate. 

In these three mechanisms, the kinetics of dissolution 
are limited by the ability of the atomic network to lo¬ 
cally rearrange in order to: (a) accommodate intermedi¬ 
ate over-coordinated atoms (Si v or O m ) or (b) enable the 
creation of pathways, or channelby enlarging preex¬ 
isting holes, thereby allowing a penetrating species (wa¬ 
ter molecule or H + ) to transfer between two sites. As 
such, the activation energy associated with such reorga¬ 
nizations takes the form of a strain energy that charac¬ 
terizes the ability of the network to locally resist, that 
is, to prevent the formation of energetic instabilities (de¬ 
parture from an equilibrium configuration) that are in¬ 
duced due to a dissolution process. Then, it becomes 
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Figure 4. Activation energy of dissolution, (a) Dissolution rates of glassy silica and o-quartz at pH 13 as a function of the 
inverse temperature 10 3 /T. The data are fitted by an Arrhenius-like law K — Ko exp(— n c Eo/RT), with Ko = 6.4 10 7 /xmol/m 2 /s 
and Ea = 52 kJ/mol for glassy silica, and Ko = 1.6 10 14 /xmol/m 2 /s and Ea = 101 kJ/mol for o-quartz (dashed lines), (b) The 
activation energy of the dissolution process for glassy silica and o-quartz at pH 13, with respect to the number of constraints 
per atom n c . The data are compared with the activation energies of sodium silicate (self-diffusiorP^ and conductiorP^ 
potassium silicate (self-diffusiorr 5 ^ and conduction 36 of K) fused silica (diffusion of water 37 ), and quartz (diffusion of w; 

The dashed line shows the effective activation energy predicted from the number of constraints per atom, following Ea = n c Eo , 
where Eq is the activation energy required to break a given constraint per atom present in the atomic network. 



clear that the development of a strain energy, and its fol¬ 
low on dissipation, is linked to the number of constraints 
per atom n c , since each constraints acts as a little spring 
between the atoms. Therefore, n c characterizes the steric 
effect imposed by the architecture of the atomic network, 
which prevents the accommodation of defects (i.e., over¬ 
coordinated atom or enlarged hole to create a channel). 
This picture is in line with result from density functional 
theory, which have shown that the activation energy of 
the hydrolysis of the inter-tetrahedra bridging oxygen in¬ 
creases with the extent of the steric constraints imposed 
by the network, that is, its connectivity and therefore, 
rigidit}E2. 


IV. CONCLUSIONS 

Our results demonstrate that the dissolution rate of 
silicate-based glassy and crystalline solids in caustic so¬ 
lutions is correlated to the number of topological con¬ 
straints per atom present in the network. Such atomic 
constraints are indicative of the effective activation en¬ 
ergy that needs to be overcome for dissolution to oc¬ 
cur. Moreover, n c characterizes the ease with which the 
atomic network can reorganize to accommodate the in¬ 
termediate (metastable) defects states that are formed 
over the course of the net dissolution process. The rate 
at which such reorganization can occur, appears to be 
the primary rate-controlling step during dissolution, that 
is, which controls its kinetics. This leads to the idea 
that dissolution kinetics, a response through to be con¬ 
trolled by surface properties, is paradoxically primarily 
controlled by the topology of the bulk architecture of 
the atomic network of the dissolving material. In addi¬ 
tion, it is remarkable that a number of processes that 


share a common thermally activated origin, e.g., dissolu¬ 
tion, diffusion and transport show similar trends in acti¬ 
vation energy with respect to the number of topological 
constraints per atom. If, indeed, it is well known that 
diffusion and conductivity are related to each other via 
the Nernst-Einstein equation, it is likely that dissolution 
shares a similar basis, although perhaps better described 
by a Butler-Volmer type of formulation^. Neverthe¬ 
less, our results suggest that dissolution, diffusion, and 
conduction could potentially be treated within the same 
atomistic, and thermodynamic framework - enabling a 
priori estimation of such phenomena. This has impor¬ 
tant implications for better understanding and modeling 
processes ranging from cement hydration, to the chemi¬ 
cal durability of glasses used to encapsulate radioactive 
wastes, to rate controls on silicate mineral dissolution 
that accompanies ocean acidification 43 . 


ACKNOWLEDGMENTS 

The authors acknowledge full financial support for this 
research provided by: The U.S. Department of Trans¬ 
portation (U.S. DOT) through the Federal Highway Ad¬ 
ministration (DTFH61-13-H-00011), the National Sci¬ 
ence Foundation (CMMI: 1066583 and CAREER Award: 
1235269), The Oak Ridge National Laboratory oper¬ 
ated for the U.S. Department of Energy by UT-Battelle 
(LDRD Award # 4000132990), and the University of Cal¬ 
ifornia, Los Angeles (UCLA). Access to computational 
resources was provisioned by the Physics of AmoRphous 
and Inorganic Solids Laboratory (PARISlab), the Labo¬ 
ratory for the Chemistry of Construction Materials (LC 2 ) 
and the Institute for Digital Research and Education 
(IDRE) at UCLA. This research was conducted in: Labo- 






6 


ratory for the Chemistry of Construction Materials (LC 2 ) 
and Physics of AmoRphous and Inorganic Solids Labo¬ 
ratory (PARISlab) at UCLA. The authors gratefully ac¬ 
knowledge the support that has made these laboratories 
and their operations possible. The contents of this paper 


reflect the views and opinions of the authors, who are 
responsible for the accuracy of the datasets presented 
herein, and do not reflect the views and/or policies of 
the funding agencies, nor do the contents constitute a 
specification, standard or regulation. 


* Contact: bauchy@ucla.edu: 

Homepage: http://mathieu.bauchy.com 

t Contact: gsant@ucla.edu; 

Homepage: http://www.lcc-ucla.edu 

1 Julian R. Jones, “Review of bioactive glass: From Hench 
to hybrids,” Acta Biomaterialia 9, 4457-4486 (2013) 

2 K.b. Harvey and C.d. Litke, “Model for Leaching Behavior 
of Aluminosilicate Glasses Developed as Matrices for Im¬ 
mobilizing High-Level Wastes,” Journal of the American 
Ceramic Society 67, 553-556 (1984) 

13 Bernd Grambow, “Nuclear Waste Glasses - How Durable?” 
Elements 2, 357-364 (2006) 

4 H. F. W. Taylor, Cement Chemistry (Thomas Telford, 
1997). 

5 Celine Cailleteau, Frederic Angeli, Frangois Devreux, 
Stephane Gin, Jacques Jestin, Patrick Jollivet, and Olivier 
Spalla, “Insight into silicate-glass corrosion mechanisms,” 
Nature Materials 7, 978-983 (2008) 

6 Chris Anbeek, “Surface roughness of minerals and implica¬ 
tions for dissolution studies,” Geochimica et Cosmochimica 
Acta 56, 1461-1469 (1992) 

‘ Susan L. Brantley, “Kinetics of Mineral Dissolution,” in 
Kinetics of Water-Rock Interaction , edited by Susan L. 
Brantley, James D. Kubicki, and Art F. White (Springer 
New York, 2008) pp. 151-210. 

8 W. A. Lanford, K. Davis, P. Lamarche, T. Laursen, 
R. Groleau, and R. H. Doremus, “Hydration of soda- 
lime glass,” Journal of Non-Crystalline Solids 33, 249-266 
(1979) 

y William H. Casey and Henry R. Westrich, “Control of 
dissolution rates of orthosilicate minerals by divalent 
metal-oxygen bonds,” Nature 355, 157-159 (1992) 

10 C. Andre Ohlin, Eric M. Villa, James R. Rustad, and 
William H. Casey, “Dissolution of insulating oxide mate¬ 
rials at the molecular scale,” Nature Materials 9, 11-19 
( 2010 ) 

11 James P. Hamilton, Carlo G. Pantano, and Su¬ 
san L. Brantley, “Dissolution of albite glass and crystal,” 
Geochimica et Cosmochimica Acta 64, 2603-2615 (2000) 

12 Roland Hellmann, Stephane Cotte, Emmanuel Cadel, 
Sairam Malladi, Lisa S. Karlsson, Sergio Lozano-Perez, 
Martiane Cabie, and Antoine Seyeux, “Nanometre-scale 
evidence for interfacial dissolution-reprecipitation control 
of silicate glass corrosion,” Nature Materials 14, 307-311 
(2015) 

id James P. Hamilton and Carlo G. Pantano, “Effects of 
glass structure on the corrosion behavior of sodium- 
aluminosilicate glasses,” Journal of Non-Crystalline Solids 
Practical Implications of Glass Structure, 222, 167-174 
(1997) 

14 James P. Hamilton, Susan L. Brantley, Carlo G. Pan¬ 
tano, Louise J. Criscenti, and James D. Kubicki, “Dis¬ 
solution of nepheline, jadeite and albite glasses: toward 
better models for aluminosilicate dissolution,” Geochimica 


et Cosmochimica Acta 65 , 3683-3702 (2001) 

10 J. C. Phillips, “Topology of covalent non-crystalline solids 
.1. Short-range order in chalcogenide alloys,” Journal of 
Non-Crystalline Solids 34, 153-181 (1979) 

1D M. F. Thorpe, “Continuous deformations in random net¬ 
works,” Journal of Non-Crystalline Solids 57 , 355-370 
(1983), 

L< J. C. Mauro, “Topological constraint theory of glass,” 
American Ceramic Society Bulletin 90 , 31-37 (2011). 

18 J. Clerk Maxwell, “L. On the calculation of the equilibrium 
and stiffness of frames,” Philosophical Magazine Series 4 
27,294-299 (1864) 

iy Aditya Kumar, Jason Reed, and Gaurav Sant, “Vertical 
Scanning Interferometry: A New Method to Measure the 
Dissolution Dynamics of Cementitious Minerals,” Journal 
of the American Ceramic Society 96 , 2766-2778 (2013) 

zu Antonio C. Lasaga and Andreas Luttge, “Variation of 
Crystal Dissolution Rate Based on a Dissolution Stepwave 
Model,” Science 291 , 2400-2404 (2001) 

21 Patricia M. Dove, Nizhou Han, and James J. De Yoreo, 
“Mechanisms of classical crystal growth theory explain 
quartz and silicate dissolution behavior,” Proceedings 
of the National Academy of Sciences 102 , 15357-15362 
(2005) 

James Patrick Hamilton, Corrosion behavior of sodium alu¬ 
minosilicate glasses and crystals (1999). 

A. K. Varshneya, Fundamentals of Inorganic Glasses (Aca¬ 
demic Press Inc, 1993). 

24 M. Micoulaut, A. Kachmar, M. Bauchy, S. Le Roux, 
C. Massobrio, and M. Boero, “Structure, topology, rings, 
and vibrational and electronic properties of Ge_{x}Se_{l- 
x} glasses across the rigidity transition: A numerical 
study,” Physical Review B 88, 054203 (2013), 

25 Mathieu Bauchy, Mohammad Javad Abdolhosseini Qomi, 
Christophe Bichara, Franz-Josef Ulm, and Roland J.- 
M. Pellenq, “Rigidity Transition in Materials: Hardness 
is Driven by Weak Atomic Constraints,” Physical Review 
Letters 114 , 125502 (2015) 

X. L. Yuan and A. Cormack, “Si-O-Si bond angle and tor¬ 
sion angle distribution in vitreous silica and sodium sili¬ 
cate glasses,” Journal of Non-Crystalline Solids 319 , 31-43 
(2003) 

M. Bauchy and M. Micoulaut, “Atomic scale foundation 
of temperature-dependent bonding constraints in network 
glasses and liquids,” Journal of Non-Crystalline Solids 357 , 
2530-2537 (2011) 

Y. Vaills, T. Qu, M. Micoulaut, F. Chaimbault, and 
P. Boolchand, “Direct evidence of rigidity loss and 
self-organization in silicate glasses,” Journal of Physics- 
Condensed Matter 17 , 4889-4896 (2005) 

zy Ye Xiang, Jincheng Du, Morten M. Smedskjaer, and 
John C. Mauro, “Structure and properties of sodium alu¬ 
minosilicate glasses from molecular dynamics simulations,” 
The Journal of Chemical Physics 139 , 044507 (2013) 



7 


30 M. Bauchy, “Structural, vibrational, and elastic properties 
of a calcium aluminosilicate glass from molecular dynamics 
simulations: The role of the potential,” The Journal of 
Chemical Physics 141 , 024507 (2014) 

M. Bauchy, “Topological constraints and rigidity of net¬ 
work glasses from molecular dynamics simulations,” Amer¬ 
ican Ceramic Society Bulletin 91 , 34-38A (2012) 

^ M. Bauchy, M. Micoulaut, M. Celino, S. Le Roux, 
M. Boero, and C. Massobrio, “Angular rigidity in tetrahe¬ 
dral network glasses with changing composition,” Physical 
Review B 84 , 054201 (2011) 

^ B. C. Bunker, “Molecular mechanisms for corrosion of sil¬ 
ica and silicate glasses,” Journal of Non-Crystalline Solids 
Proceedings of the First PAC RIM Meeting on Glass and 
Optical Materials, 179 , 300-308 (1994) 

35 Siting Zhang and Yun Liu, “Molecular-level mechanisms 
of quartz dissolution under neutral and alkaline conditions 
in the presence of electrolytes,” Geochemical Journal 48, 
189-205 (2014) 

130 Gunther Heinz Frischat, Ionic diffusion in oxide glasses 
(Trans Tech Publications, 1975). 

36 K. L. Ngai and S. W. Martin, “Correlation between the ac¬ 
tivation enthalpy and Kohlrausch exponent for ionic con¬ 
ductivity in oxide glasses,” Physical Review B 40 , 10550- 
10556 (1989) 

A. J. Moulson and J. P. Roberts, “Water in silica glass,” 
Trans. Faraday Soc. 57 , 1208-1216 (1961). 

38 E. W. Shaffer, J. Shi-Lan Sang, A. R. Cooper, and A. H. 
Heuer, “Diffusion of tritiated water in a-quartz,” Sam by 
BJ Gjdleto H, A Yoder, RA 1, 131 (1974). 


39 Bruno MJ Smets and T. P. Lommen, “Role of molecular 
water in the leaching of glass,” Physics and Chemistry of 
Glasses 24 , 35-6 (1983). 

40 M. Bauchy and M. Micoulaut, “From pockets to channels: 
Density-controlled diffusion in sodium silicates,” Physical 
Review B 83, 184118 (2011) 

41 Alexander Pelmenschikov, Helene Strandh, Lars G. M. 
Pettersson, and Jerzy Leszczynski, “Lattice Resistance to 
Hydrolysis of Si-O-Si Bonds of Silicate Minerals: Ab Ini¬ 
tio Calculations of a Single Water Attack onto the (001) 
and (111) /3-Cristobalite Surfaces,” The Journal of Physical 
Chemistry B 104 , 5779-5783 (2000) 

Kai Kristiansen, Markus Valtiner, George W. Greene, 
James R. Boles, and Jacob N. Israelachvili, “Pressure 
solution - The importance of the electrochemical surface 
potentials,” Geochimica et Cosmochimica Acta 75 , 6882- 
6892 (2011) 

James C. Orr, Victoria J. Fabry, Olivier Aumont, Lau¬ 
rent Bopp, Scott C. Doney, Richard A. Feely, Anand 
Gnanadesikan, Nicolas Gruber, Akio Ishida, Fortunat 
Joos, Robert M. Key, Keith Lindsay, Ernst Maier- 
Reimer, Richard Mat ear, Patrick Monfray, Anne Mouchet, 
Raymond G. Najjar, Gian-Kasper Plattner, Keith B. 
Rodgers, Christopher L. Sabine, Jorge L. Sarmiento, 
Reiner Schlitzer, Richard D. Slater, Ian J. Totterdell, 
Marie-France Weirig, Yasuhiro Yamanaka, and Andrew 
Yool, “Anthropogenic ocean acidification over the twenty- 
first century and its impact on calcifying organisms,” Na¬ 
ture 437 , 681-686 (2005) 


